Particle tunneling through a polarizable insulator 
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The tunneling probability between two leads connected by a molecule, a chain, a film, or a bulk 
polarizable insulator is investigated within a model of an electron tunneling from lead A to a state 
higher in energy, describing the barrier, and from there to lead B. To describe the possibility of 
energy exchange with excitations of the molecule or the insulator we couple the intermediate state 
to a single oscillator or to a spectrum of these, respectively. In the single-oscillator case we find for 
weak coupling that the tunneling is weakly suppressed by a Debye- Waller-type factor. For stronger 
coupling the oscillator gets stiff and we observe a suppression of tunneling since the effective barrier 
is increased. The probability for the electron to excite the oscillator increases with the coupling. In 
the case of a film, or a bulk barrier the behavior is qualitatively the same as in the single oscillator 
case. An insulating chain, as opposed to a film or a bulk connecting the two leads, shows an 
orthogonality catastrophe similar to that of an electronic transition in a Fermi gas. 

PACS numbers: 73.63.-b, 73.23.Hk, 72.20.Dp 



There is a vast interest in the electronic properties of 
molecular systems or tunnel junctions due to their po- 
tential use for electronic applications. Thereby, a tun- 
nel junction might be a microscopic or macroscopic one-, 
two- or three dimensional insulator connecting two metal 
leads. As long as the molecule or the tunnel junction is 
not excited by the electron the conductance can be phe- 
nomenologically described by introducing a matrix ele- 
ment for tunneling through the barrier. However, the 
possibility to excite the barrier opens new paths for the 
electron, thus, enhancing the conductance. In recent 
experiment on molecular junctionil'EI the conductance 
shows several steps with increasing bias voltage between 
the two leads underlining the inmortance of vibrations 
excited by the tunneling clcctrontl. We are solely inter- 
ested in the influence of vibrations on the conductance 
and we do not restrict ourselves to weak coupling between 
the vibrations and the tunneling electron. In detail we in- 
vestigate the probability for a tunneling electron to excite 
the molecule/tunnel junction. Therefore we introduce an 
intermediate state in the junction which itself is coupled 
to either a single oscillator to mimic the excitations of 
a molecule or to a spectrum of oscillators to mimic pho- 
tonic excitations of a tunnel junction. A similar approach 
allowing two paths to a inal state for a tunneling parti- 
cle was used by HarrisoriQ to investigate how the coupling 
destroys the coherence between two paths. 



A Single Oscillator — The electron starts in lead A 
from an initial state |0) with energy Eq. This state is 
coupled to the intermediate electron state |1) with en- 
ergy El > Eq by the tunneling matrix element Vqi- The 
intermediate state is also coupled by a matrix element 
Via/VM to M final states \a) with a > 2 in lead B. For 
sufficiently small tunneling matrix elements the transi- 
tion probability Po,f for the electron tunneling from |0) 
through |1) to a final state \a) can be obtained by simple 



second-order time-dependent perturbation theorytl as 
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with e = El — Eq ^ Vbi,ViQ. For the summation over 
all final electron states we assumed a uniform density of 
states of the M states within a band of bandwidth h in 
lead B. 

An interaction Hi between the intermediate state and 
the excitations of a molecule can be mimicked by cou- 
pling the state |1) to a single harmonic oscillator: Hi = 
7|l)(l|i: where 7 is a coupling constant and x is the dis- 
placement coordinate of the oscillator. The quantum me- 
chanical oscillator states when the electron is in one of 
the leads are denoted as \4>i)^ where i is the excitation 
number, whereas is the oscillator state with i excita- 
tions when the electron is in the intermediate state. The 
total state of the problem can be expressed by the states 
|0, and |2,(/)j) with energy iJj^^o.i.a + (the 

zero-point energy is absorbed in the En). 

We consider only initial states without excitations. 
The probability Poo.fn for tunneling from state |0) with 
zero excitations to a final state with n excitations in the 
oscillator is modified from Eq. (nl) as 



-Poo,fn — -Fb, 
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where we sum over all possible numbers i of virtual exci- 
tations in the oscillator in the intermediate state. 

The oscillator states tpi^x) = (j)i{x ~ xq) are shifted by 
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since an electron in the intermediate state shifts the po- 
tential minimum position of the oscillator. With the zero 
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point fluctuation 
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of the harmonic osciUator we obtain for the overlaps 



= exp 



E 



4ag 



(5) 



T2T Xo 



- n + j)\{n - i)\ Vflo 



Assuming that the osciUator is in the ground state before 
and after tunnehng the expressions simphfy substantiaUy 
and we obtain for the tunnehng probabihty 
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exp(-2W^) , (6) 
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The sum runs over the number i of virtual excitations 
in the intermediate state. The contributions peak for 
i W virtual excitations. In the weak-coupling limit, 
1, tunneling without virtual excitations, i = 0, 
dominates and the tunneling is influenced by a Debye- 
Waller type term ~ exp(— 2VF). For strong coupling, 
3> 1, the contribution to the sum peaks a,i i ^ W and 
noting /il ~ we obtain 
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which numerically turns out to be very precise when 
W ^ 1. The tunneling electron leaves the oscillator 
unmoved but its potential minimum position is shifted 
putting the oscillator in an excited state with an excita- 
tion energy imw^Xg = WTilo. Accordingly, the barrier is 
effectively raised by this energy. 

Naturally, we would assume that the crossover at 
W ^ 1 coincides with the region where Poo.fo falls from its 
value without coupling to an oscillator to a strongly sup- 
pressed probability. However, Eq. (||) shows that even for 
14^ ^ 1 the coupling to the oscillator only influences the 
tunneling probability substantially for Whuj > e. Our 
perturbative approach limits us to energies smaller than 
the barrier heights. Therefore we cannot investigate the 
strong coupling limit Wfioj > e and we are confined to 
oscillator frequencies u < e/h. 

Defining a number rir with 



rhuj — 



(9) 



we find that Poo.fo is a function of W/uj. instead of W 
and we are restricted to the effective weak coupling limit 

W/Ur < 1. 



If we interpret the virtual excitation energy Whco as 
the energy put into the oscillator during tunneling we 
expect that in the resonant case when W = 1 the prob- 
ability for exciting the oscillator will be maximal. How- 
ever, we saw that for huj < e the influence on the electron 
from the coupling to the oscillator increases with W/ur- 
We numerically find that the probability Poo, 21 for excit- 
ing the oscillator with a single excitation increases lin- 
ear with W/rir (quadratic with the coupling 7 between 
electron and oscillator) and the probability Poo,2n de- 
creases for a higher number of excitations n. If we ex- 
pand e/(e -|- ihuj) « 1 — ihuj/e in Eq. (^) up to first order 
in fiuj/e we obtain 
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which shows the quadratic behavior with the coupling 7. 
Summing over all possible numbers n of excitations we 
get for the probability to excite the oscillator Poo.fn>i ~ 
(7ao/e)^Po,f . This approximation has the same limit as 
discussed before namely that the virtual energy Whuj of 
the mainly contributing tunneling path is smaller than 
the barrier heights e > WTiuj. Therefore this result holds 
only for couplings where (jao/e)'^ < (27itj/e)Po,f . 
Multiple Oscillators — In the next step we introduce 
a spectrum of oscillators coupled to the tunneling elec- 
tron to investigate the tunneling of an electron through 
a tunnel junction, for example a linear chain of atoms as 
illustrated in Fig. |l|. Assuming that all bonding states 
between the atoms of the chain are occupied by electrons, 
the tunneling electron hops to an anti-bonding state be- 
tween two atoms. The electron repels with force F the 
nearest neighbors resulting in a coupling 

Hi = -|1)(1| -P - (fc™+i - fc™) 

of an electron in state |1) to the displacement Sx„i-i-i and 
Sxm of the nearest neighbors. For a chain of N atoms 
and lattice constant d we obtain, using fixed boundary 
conditions at position r = and r = {N + l)d, N phonon 
modes with 6xm — y^T^Tlk ""fc sin(fcxm) and 

TL TT 

k = kn = with n=l..N . (11) 

It is convenient to take N even and the anti-bond m — 
N/2. Then the coupling between the electron and the 
phonons turns out to be 

i// = |l)(l|-^Afe„Ufc„ 



with the coupling constants 
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: n even 



In the following we use the indices k and n inter- 
changeably. The coupling constant A„ between the elec- 
tron and the mode with wave vector fc„ is proportional 
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FIG. 1: A linear chain of atoms is illustrated in which a tun- 
neling electron might hop to an anti-bond between atom m 
and m + 1 repelling both neighbors. 



to 1/VN, so in the thermodynamic limit {N — > oo) the 
coupling effects estimated by ^^^Y finite. 

In the case of a single oscillator the behavior of the 
system was determined hy W = ^(xo/ao)^. Since A„ ~ 
in the case of many phonons each single mode is 
weakly coupled to the electron for large N. However, 
the influence to the tunneling electron is now determined 
by the sum over all modes W = 5 with 
Xn — ^kn/i'muj'^) and a„ = ^h/ (ma;„). The equations 
given above are for a one-dimensional chain but they are 
easy to generalize to two and three dimensions. It is a lit- 
tle intricate to work out the coupling but essentially the 
coupling constants look the same and only the den- 
sity of phonon modes g{k) depends on the dimensionality, 
g{k) (X k'^~^. This turns out to give a striking difference 
between the one dimensional and the other cases. 

With no phonon in the initial state {ruk = 0} we find 
for the transition probability to the final state with Ik 
excitations in mode fc„ 



Pa{nik=o}.i{h} - Po,{ ■ 



(13) 



E 



where (k) and tpi^ (fc) are the wavefunctions of the un- 
shifted and shifted mode k with I excitations respectively. 
The sum sums over all possible total transient 

states with virtual excitation numbers ik for each mode. 
For each of these states the product Ylk ^'^^^ the various 
overlaps calculates the total overlap between the initial 
and the transient state times the overlap of the transient 
and the final state. 

Let us focus first on the probability Po{mk=o},i{ik=o} 
for tunneling without final state excitations. The con- 
tribution ^0 from the tunneling 'path' without virtu- 
ally exciting any phonon is given by the overlap between 
the ground states of the unshifted and shifted oscillator 
modes 
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Einstein Approximation - An Einstein approximation 
leaves us essentially with N modes with the same fre- 
quency uje and coupling constant 7. We find 
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which is identical to the single oscillator result. For the 
contributions A^ from tunneling paths involving m vir- 
tual excitations we obtain 
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The factorial 1/m! results from the fact that, once we 
chose m oscillators to put the m virtual excitation in, 
there are still m! ways of distributing the virtual excita- 
tion to the oscillators. Physically, however all these m! 
ways are identical. Technically it is convenient to sum 
over all these and we obtain mlAm which we then divide 
by the factorial to get Eq. (|l^). Note also that Am in- 
cludes only tunneling path where all m excitations are 
in m different modes. Tunneling paths with more than 
one excitation in a mode are suppressed for large N by 
factors 1/A^. The probability /o{mt=o},f{Zfc=o} for tunnel- 
ing without final state excitations thus is equivalent to 
Eq. (^ . An Einstein spectrum yields the same tunneling 
probability as a single oscillator if no mode is excited. 

The probability for tunneling through the chain and 
exciting mode j with a single excitation is given by 
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Pqq 20 is the probability for the electron tunneling 
through the chain without final state excitations calcu- 
lated without including mode j. Since a single mode does 
not change results appreciably we take -P00.20 — -Poo, 20 • 
We sum j over all modes. The exponential is unity for 
large N and we use Eq. (|lj) to find for the probability 
to excite n modes 



-Fbo.fri — —,W^ ^00, 20 
n! 



(17) 



Note that contributions from tunneling path with more 
than one excitation in a single mode are suppressed. The 
factorial n\ has the same origin as m! in Eq. (p^. For 
W < 1 tunneling through the chain without exciting it 
and without using virtual excitations dominates. How- 
ever, for W > 1 the tunneling electron uses predom- 
inantly W virtual excitations to tunnel and the chain 
has most likely W excitations once the electron tunneled 
through. This shows a difference from the single oscilla- 
tor case where in the strongly coupled regime, W > 1, 
the oscillator was not excited by the tunneling electron. 
Only for W/rir near one excitations are likely. To un- 
derstand this we note that each single mode in the many 
oscillator case is weakly coupled. The probability to ex- 
cite some is only finite because there are many modes. 
Debye Approximation - Within a Debye approximation 
we assume a spectrum of modes with linear dispersion 
uo — vk with the speed of sound v. The Debye cut-off 
ku — {N/N + l)(7r/c?) is chosen to give the correct num- 
ber of modes. It leads qualitatively to the same results as 
an Einstein approximation for a two and a three dimen- 
sional tunnel junction since their phonon spectra favor 
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high frequency modes and is replaced by an average 
frequency. We therefore focus on tunnehng through a 
chain of atoms. We further restrict the following con- 
siderations to the probability for tunneling without final 
state excitations since the probability for tunneling with 
final state excitations follows then from Eq. (|l7|). 

The chain has a constant density of states and we find 

with fcc = (Fd)'^ /{hmv^). For large N tunneling through 
the chain without the help of virtual excitations is 
suppressed by exp(— 2Vl^) ~ exp(— 2(fcc/^z?) In(Af)) = 
j^-(2k^/ko) ^ In the thermodynamic limit, N — > oo, this 
approaches zero. This behavior is similar ta.the 'orthog- 
onality catastrophe' described by AndersonQ for an elec- 
tron tunneling into a Fermi gas which has a constant 
density of states near the Fermi surface. 

We now investigate the contribution Am from a tunnel- 
ing path involving m virtual excitations. The key ques- 
tion is whether there is a dominant tunneling path (a m 
with Am maximal) and what would be its effective vir- 
tual energy J^k^k^'^k (compare Eq. ([l3|)). We suggest 
that the contribution of mode j to the effective virtual 
energy is determined by a factor (x^/a^) since the same 
factor determined in the case of tunneling with one final 
state excitation the contribution of the tunneling path 
with mode j finally excited (compare Eq. (p^)). Since 
Am involves m virtual excitations we introduce further a 
'normalizing' factor m/W and find 

For the energy of the m virtual excitations we find 

Emx'^ m 
= — hujc . (20) 

Numerically we find that this approximation holds when 
the energy of the m virtual excitations does not exceed 
the barrier heights. For large finite N we find that W 
tends logarithmically to infinity. However, for W ^ 1 
the probability for tunneling without final state excita- 
tions is simply given by the tunneling probability through 



an effective barrier (compare Eq. (||)) enhanced by the 
energy of the virtual excitations of the dominating tun- 
neling path. The dominating path involves W virtual 
excitations leading to an energy fujc of the virtual exci- 
tations which is independent of W . Thus the probability 
for tunneling through the chain without final state exci- 
tations is finite and independent of the length N . 

Note that the maximal virtual energy WTiud possi- 
bly used in a tunneling event with W virtual excitations 
increases logarithmically with the number N of modes. 
Once this energy WTilod exceeds the barrier heights e the 
numerical investigation showed that our approximation, 
Eq. (|l^), fails. However at the same point our approach 
using second order perturbation theory is questionable. 
We therefore believe that our approximative result holds 
in the thermodynamic limit and that tunneling through 
a chain without final state excitations is possible. 

In conclusion, we first investigated tunneling through 
a barrier with a coupling to a single oscillator. The re- 
sult can be characterized by a coupling strength W given 
by the square of the ratio of the classical shift of the os- 
cillator to its zero-point fluctuation. For weak coupling, 
W < 1, tunneling without final state excitations is sup- 
pressed by a factor e~^^ « 1 — 2W . For strong coupling, 
W > 1, tunneling is suppressed by an increase in the 
tunneling barrier height by an energy W times the vibra- 
tional energy hu. The probability for tunneling with final 
state excitations increases in proportion to W . We then 
introduced a barrier with coupling to a spectrum of N 
modes, with couplings proportional to 1/ \/N . It is shown 
that in an Einstein approximation, where all modes have 
the same frequency uje, the results are essentially the 
same. The same holds true for two or three dimensional 
Debye spectra where the frequency is proportional to the 
wave number. However, a one-dimensional chain shows 
an orthogonality catastrophe similar to that of an elec- 
tronic transition in a Fermi gas since W oc ln(A'^). The 
probability for tunneling without final state excitations 
and without the help of virtual excitations is suppressed 
by a factor 1/iV raised to a finite power. However, our 
investigation suggests that the probability for tunneling 
without final state excitations but with the help of virtual 
excitations stays finite for N ^ oo. 
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